Draft version November 3, 2010 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



A MODEL FOR THERMAL PHASE VARIATIONS OF CIRCULAR AND ECCENTRIC EXOPLANETS 

Nicolas B. Cowan 1 ' 2 , Eric Agol 2 , 

Draft version November 3, 2010 

ABSTRACT 

We present a semi-analytic model atmosphere for close-in exoplanets that captures the essential 
physics of phase curves: orbital and viewing geometry, advection, and re-radiation. We calibrate the 
model with the well-characterized transiting planet, HD 189733b, then compute light curves for seven 
of the most eccentric transiting planets: Gl 436b, HAT-P-2b, HAT-P-llb, HD 17156b, HD 80606b, 
WASP-17b, XO-3b. We present phase variations for a variety of different radiative times and wind 
speeds. In the limit of instant re-radiation, the light curve morphology is entirely dictated by the 
planet's eccentricity and argument of pericenter: the light curve maximum leads or trails the eclipse 
depending on whether the planet is receding from or approaching the star at superior conjunction, 
respectively. For a planet with non-zero radiative timescales, the phase peak occurs early for super- 
rotating winds, and late for sub-rotating winds. We find that for a circular orbit, the timing of the 
phase variation maximum with respect to superior conjunction indicates the direction of the dominant 
winds, but cannot break the degeneracy between wind speed and radiative time. For circular planets 
the phase minimum occurs half an orbit away from the phase maximum — despite the fact that the 
coolest longitudes are always near the dawn terminator — and therefore does not convey any additional 
information. In general, increasing the advective frequency or the radiative time has the effect of 
reducing the peak-to-trough amplitude of phase variations, but there are interesting exceptions to 
these trends. Lastly, eccentric planets with orbital periods significantly longer than their radiative 
time exhibit "ringing" whereby the hot spot generated at periastron rotates in and out of view. The 
existence of ringing makes it possible to directly measure the wind speed (the frequency of the ringing) 
and the radiative time constant (the damping of the ringing). 
Subject headings: methods: data analysis — (stars:) planetary systems — 



1. INTRODUCTION 

Efforts to model the atmospheres of exoplanets fall 
broadly into two categories: 1) radiative transfer models 
that make predictions about the planet's spectral charac- 
teristics after making assumptions about the bulk energy 
trans port in the atmosphere (numerical simulations in- 
clude llro et al.l 120051: ISeager et al.l 120051: IBarman et al.l 
120051: iBurrows et"aLl l200a iFortnev et al.l 120081 while 



Hansen 2008 developed a series of analytic radiative 
transfer models, and Madhusudhan & Seager 2009 
introduced a Monte Carlo approach for fitting models to 
observed spectra); 2) hydrodynamic models that make 
predictions about energy recirculation given assumptions 
about radiative transfer. The hydrodynamic models 
run with a relatively simple treatment of radiative 
transfer may then be p ost-processed w i th mo re detailed 
radiative transfer (e.g.. [Burrows et al.ll2010l) . Roughly 
in increasing order of sophistication, the hydrodynamic 
models consist of solving the: sh allow water equations 
(jLangton fc La ughlin 2007, 2008 ), equ iv alent barotropic 
equat ions ( Cho et all 12003). | 2008at [Rauscher et alJ 
20081). primitive equ a tions (IShowman fc Guillotl 12002 



Cooper fc Showman 120051 : iShowman et al 
Menou fc Rausrlu^ [200 



iRauscher fc Menoul 
.Thrastarson fc Cho I 12010ft . Euler's e quations 
(Burkertetal .1 120051 : IDobbs-Dixon fc Lml [2008) , or the 
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full Navier Stokes equation ( Dobbs-Dixon et ail 120101) 
Recen tly, IShowman et al.l ( 2009a|) and iLewis et al.l 
(2010) have combined approaches 1 and 2 by simultane- 
ously solving the radiative transfer and hydrodynamical 
equations. With the exceptions o f the h yd rodynamical 
model s of lLangton fc Laughlin! (|2008| ). ILewis et al.1 
(2010) and iBeaulieu et al.1 (|2010|). as well as the radia- 
tive transfer models of llro fc Deming (2010), however, 
these theoretical efforts have been directed towards 
understanding exoplanets on circular orbits, for which 
the planet's power budget is constant. 

Despite the difficulties in numerically modeling exo- 
planet atmospheres, the physical parameters that govern 
their large-scale behavior are thought to be: the radiative 
timescale (how quickly a parcel of gas radiates away its 
thermal energy), and the advective timescale (how long 
it takes for a parcel of gas to move from the planet's day- 
side to its night-side). In detail, not every joule of energy 
is absorbed at the same depth (let alone re-radiated from 
that depth) and different parcels of gas will take differ- 
ent amounts of time to move from the day-side to the 
night-side. Nevertheless, lacking detailed observations 
of these planets, such a simple approach complements 
more detailed simulations by indicating — at very little 
computational cost — which ingredients are necessary to 
properly explain observed behavior. 

The radiative timescale is primarily determined by 
the depths at which incident energy from the host star 
is being absorbed and emitted (which are not neces- 
sarily the same). Radiative transfer models suggest 
that this occurs at pressures of 10 _3 -10° bars, where 
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the radiative time scales are roughly 10 4 -10 5 seconds 
(ISeaeer et all 120051; jIro~eTaII [200a iFortnev et all [20051: 
iBarman et al.l 120051 : iBurrows et al. 2006). Advection 
timescales are a function of planetary size and wind 
speed. Different hydrodynamic simulations put the 
wind speeds at roughly a tenth of the sound speed 
at the slow end, and transonic at the high end: ~ 
10 2 -10 3 m/s (IShowman fe Guillotll200l ICho et al.l | 2003 : 



Burkert et al.l 120051: ICooper fc Showman! I2005L 120061 : 
Langton fc Laughlinl 120071: iThrastarson fc Choi I20T0I ) 
Such wind speeds lead to advective times of ~ 10 5 s for a 
Jupiter-sized planefl Since the radiative and advective 
timescales are of the same magnitude, it is not immedi- 
ately clear how efficient energy recirculation should be 
for short-period exoplanets. Indeed, hydrodynamic sim- 
ulations performed by different groups (see above) do not 
always agree. 

Fortunately, there is a growing body of observa- 
tions that constrain theoretical models for short-period 
planets. Observations of secondary eclipses in ex- 
oplanetary systems have made it possible to mea- 
sure the integ rated day-side flux of short-period ex- 
oplan ets (e.g., iCharbonneau et ail 120051 : iDeming et al.l 
2005). Observations of phase variations have made 
it possible to measure the brig htness of plan- 
ets a s a function of l ocal time ([Harrington et all 
2001 ICowan et al l [20071: IKnutson et all 1200112 00%: 
Snellen et al.l 120091: iBqrucki et all 120091: knutson et all 
2009al: iCrossfield et all 12010ft . Ideally, there would be 



analytic models to bridge the gap between the detailed 
simulations and these cutting-edge phase observations. 
We presented a summary and analysis o f the photomet- 
ric obs ervations of 24 transiting planets in lCowan fc Ago] 
(2010). We concluded that short-period planets exhibit 
a wide range of recirculation efficiencies but are all con- 
sistent with low Bond albedos (^4 < 0.35). The model 
used in that paper, however, was too simple to predict 
the shape of phase variations , only their amplitude . 

Phase function mapping (|Cowan fc Agoll 12008) is a 
largely model-independent technique, but it rests on 
the assumption that large-scale weather patterns on the 
planet do not shift with respect to the sub-stellar point. 
There is theoretical evidence that hot Jupiter atmo- 
spheres reach a steady-state diurnal heating pattern over 
a wide range of pressures that incl ude the likely mid-IR 
photo sphere (2.5 mbar - 19.6 bar, ICooper fc Showman! 
1200511 But this lac k of variability may be an artifact 
of initial conditions (Thra starson fc C ho 2010]). Never- 
theless, repeat observations of HD 189733b with Spitzer 
indicate that the 8 micron day-side variability is less 
than 2.7% and night-side variabil ity is less than 17% 
(|Agol et all 12001 1 Agol et al.ll2010h . Thus, the steady- 
state hypothesis seems to be in good shape for the large- 
scale climate of planets on circular orbits, but it is im- 
possible for this condition to hold true on planets with 
eccentric orbittQ. 

Here wc develop an analytic energy transport model — 
complementary to hydrodynamic models — which offers 

3 Unlike the radiative timescale, which depends sensitively on 
depth, the wind velo city decreases by only a fa ctor of 3 between 
2.5 mbar and 20 bar (Cooper & Showman 2005). 

4 Unless, of course, the radiative timescale is longer than an 
orbital period, in which case the planet would not exhibit any 
phase variations (sec Equation l9t . 



important insight into the bulk energy flow in the atmo- 
spheres of hot Jupiters, and its impact on the primary 
observable: thermal phase variations. In § 2 we intro- 
duce our analytic model and the associated differential 
equation (DE). We describe general solutions and useful 
limits of the DE in § 3. In § 4 we explain how we com- 
pute disc-integrated light curves. We compute sample 
light curves in § 5, discuss our model assumptions and 
how they likely impact our results in § 6, and state our 
conclusions in § 7. 

2. ENERGY TRANSPORT MODEL 

In this section, we describe a one-dimensional (longi- 
tudinal) analytic model of the heating and cooling of a 
planet. For now we simply state our simplifying assump- 
tions and develop the equations governing our model. In 
§ 6 we discuss these assumptions explicitly. We assume 
that all of the energy is reflected, absorbed and emitted 
by a single layer, which should be taken to approximate 
the planet's photosphere, or emitting surface. 

The temperature and — for planets in edge-on orbits — 
the visibility of parcels of gas is greatest near the equa- 
tor. Since we are interested only in disc-integrated light 
curves, we are principally concerned with the motion 
of the atmosphere near the equator. We assume that 
the equatorial jet stream has a constant angula r velocity 
with latitude (i.e. solid body rotation. See also llro et all 
|2005j). We define this constant effective rotation period 
for the planet, aj rot , in an inertial frame. Note that 
from the perspective of a parcel of gas, it doesn't mat- 
ter whether this motion is a result of actual rotation of 
the planet, winds, or both. Winds are usually defined 
with respect to a rotating core. It is difficult (currently 
impossible) to determine the rotation rate of a gaseous 
exoplanet's core, uj cotc . Observationally, the only acces- 
sible quantity is the sum of the planet's internal rotation 
and the winds, which we refer to as the planet's effective 
rotation: Lu rot — cu COIC + w W md- Note that the effective 
rotation of the planet can be trivially converted to an 
effective rotational period, P ro t = 27r/u; ro t, or an equato- 



rial parcel velocity in an inertial frame, v g 



R v U!r. 



For the energy budget of the parcel of gas, what mat- 
ters is not the rotation frequency but rather the related 
advective frequency, w a dv(i) = ^rot — <-^oi-b(t), where w or b 
is the angular velocity of the planet about its host star, 
which varies with time for an eccentric orbit. The ad- 
vective frequency is a measure of how often the parcel 
of gas is subjected to stellar forcing. A planet with a 
tidally locked atmospher^ has u; a d v = 0; super-rotating 
winds (in the same sense as the planet's orbital mo- 
tion) have w a( jv > 0; sub-rotating winds have w a d v < 0. 
The orbital angular velocity is a maximum at perias- 
tron: w max = (2 7 r/P)(l-e)- 3 / 2 (l-(-e) 1 / 2 , where P is the 
planet's orbital period. Note that w max is not the same 
as the pseudo-synchronous rotational frequency (see § 6 
for discussion on this). In this study, we will often ex- 
press Lu rot in units of w max since this immediately indi- 
cates whether the planet's atmosphere is super-rotating 
or trailing when its day-side is receiving the largest inci- 
dent flux. 

2.1. Energy Transport Differential Equation 
5 As opposed to a tidally locked core. 
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The flux incident on a planet changes as a function of 
time if the planet's orbit is not circular: 



F(t) 



W) 



(1) 



where a is the Stefan-Boltzmann constant, T c g and i?* 
are the star's effective temperature and radius, and r(t) 
is the planet-star distance. 

We assume that parcels of gas move only in the E- 
W direction, not latitudinally, and we adopt the La- 
grangian approach of tracking an individual parcel of gas 
and its temperature. The flux absorbed by the parcel is 
(1 — A)F(t) sin#max(cos 0), where A is the planet's 
Bond albedo (which we assume to be constant), 9 is the 
co-latitude of the parcel (9 = at the north pole, ir/2 
at the equator, and ir at the south pole), and $ is the 
angle from the sub-stellar longitude (a.k.a. local time: 
<f> = — 7r/2 at the dawn terminator; $ = at the substel- 
lar meridian; $ = ir/2 at the dusk terminator; <I> = ir at 
the antistellar meridian). 

We further assume that all of the heat absorbed by the 
parcel goes into increasing its temperature, rather than 
doing mechanical work. The warmed gas then radiates 
as an ideal blackbody with temperature T (This is sim- 
ply Newtonian cooling, a scheme used in many of the 
hydrodynamic simulations listed in § 1). The differential 
equation governing the temperature of this parcel of gas 
as a function of time is: 

— = — ((1- A)F(t) sin0max(cos$(i),O) - crT 4 ) , 
at Cfi 

(2) 

where c/j = pCpH is constant everywhere on the planet, 
provided that the mass-density, p, specific heat capacity, 
Cp, and the thickness of the parcel, H are planet-wide 
constants. Since all the incident flux is absorbed (and 
eventually re-radiated) in a single layer, one must take 
this layer to be many atmospheric scale-heights thick: 
including the optical and infrared photospheres. Note 
that even if the dynamical layer of the atmosphere is a 
dozen scale heights thick, it is still much smaller than the 
planetary radius. 

Setting the left side of Equation [2] to zero, we define 
the equilibrium temperature for a parcel of gas: 



(1 - A)F(t) sin6>max(cos$(£),0) 



(3) 

A parcel cools if its temperature is greater than T eq ; it 
heats up if its temperature is below T cq . The temper- 
ature extrema of the parcel (and any inflection points) 
thus occur when its temperature reaches T cq . Note that 
latitude enters the differential equation in the same way 
as Bond albedo. This means that locating a parcel far- 
ther from the equator is precisely equivalent to increasing 
its Bond albedo: (1 — A) «-> sin0. 

The energy budget of an eccentric planet is dictated 
by what happens near periastron. Noting that the time 
variability of T cq is contained in $(i) and F(t), we define 
the fiducial tempera ture along the sub-stellar meridian a t 
periastron (see also lHan sen 20081: iCowan fc Agoll201dO : 



T = T eS (l-A) 1 ' A sm 



V4 Q ; n l/4, 



R* 



a(l-e) 



(4) 



Following EoeLalJ (|2005t) and iSeager et all (j2005ft . we 
define the radiative time constant in the vicinity of the 
photosphere as 

c h loTl (5) 



7"rad 



We define the dimensionless temperature, T = T/Tq, 
and dimensionless tim^l t = t/r ra( j, and rewrite Equa- 
tion [2] as 

where I = max(cos 0)(a(l — e)/r(t)) 2 is the normal- 
ized intensity of radiation experienced by the parcel of 
gas. 

3. ENERGY TRANSPORT SOLUTIONS 

Equation[6]is a first-order, non-linear ordinary differen- 
tial equation (D.E.) and does not have an analytic solu- 
tion on the planet's day-side. On the planet's night-side, 
1 = 0, and the equation admits an analytic solution: 



1/3 



T = T, 



dusk 



3(i-t dusk )T d 3 usk + l 



(7) 



where Xd us k and idusk are the dimensionless temperature 
and time at sunset. 

We can — and eventually do — solve the D.E. numer- 
ically, but in the remainder of this section we consider 
simplifying limits that offer insight into the model be- 
havior. Readers primarily interested in the calculated 
light curves may skip ahead to § 4. 

3.1. The Limits of Low and High r rac i 

The limiting case of r ra( j = corresponds to a planet 
where all the flux is re-radiated from the day-side, re- 
gardless of how rapidly the gas parcels are moving. Ev- 
ery parcel of gas in such a case is always right at its 
equilibrium temperature: T($,6>,t) = T oq (<£>, i), given 
by Equation 3. 

Given our model assumptions, the r ra d 3> P mt limit 
corresponds to a planet that emits radiation uniformly 
at all longitudes. For a planet on an eccentric orbit, 
the limit comes in two varieties: if P ro t <C T ra( j <C P, 
the planet instantly re-radiates all incident energy, but 
uniformly in longitude. This may be thought of as a 
planet with strong zonal winds, leading to efficient day- 
night heat circulation: 



r '"'" ° T "[W) 



1/2 



1/4 



sin 1 / 4 9(1 - A) i/4 . (8) 



- 4^/4 



If P ro t <F< r ra d , the radiative time is so long that one 
probably cannot neglect latitudinal (meridional) trans- 
port; it is expedient to assume perfect latitudinal trans- 
port. The planet will emit uniformly everywhere on the 
planet and at the same rate all the time: 



T = T, 



cir 



2a 



1/2 



(l-A) 1 / 4 (l-e 2 ) 



-1/8 



(9) 



6 In the T ra< j — y limit, the dimensionless time is p oorly defined, 
but Equation [2] can be trivially solved, as shown in ij l3.ll 
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3.2. Circular Orbit Limit 

In the special case of circular orbits (the majority of 
short-period exoplanets, for good reason), there are many 
useful analytic approximations that we can make to sim- 
plify the problem. Furthermore, the simpler circular case 
offers intuition into the behavior of our model. In the cir- 
cular limit, Equation [6] can be rewritten as: 

^ = - (max(cos $, 0) - f 4 ) , (10) 

where e = T ra d^adv is a dimensionless constant quan- 
tifying the planet's energy recirculation efHciencjQ, and 

A day in the life of a parcel of gas proceeds as shown 
in Figure [TJ Note that we have included the sin 1 / 4 6 
factor in the y-axis, so that we may plot on the same fig- 
ure the heating curves for parcels at different latitudes. 
The two effects of a high e are 1) a delay in the time 
of maximum temperature, and 2) higher night-time tem- 
perature. Although the observed light curve depends on 
both the radiative and advective timescales, the heating 
pattern of a parcel of gas depends only on their ratio, 
e. The night-time temperature is largely independent of 
latitude, 9, but depends sensitively on e. The maximum 
temperature reached by a parcel, on the other hand, de- 
pends sensitively on its latitude but only weakly on e. 
Because of the equivalency of latitude and albedo, the 
the dashed lines in Figure Q] can be thought of as the 
heating curves for equatorial parcels of gas, but with 
A w 30%: albedo has a more important impact on the 
day-side heating pattern than on the night-side cooling. 
Finally, the delay between a parcel passing through the 
sub-stellar longitude and reaching its maximum temper- 
ature, $ m ax, depends more sensitively on e than does 
the maximum temperature reached, T max . In the e — > oo 
limit, Tsin 1//4 # = (I/71") 1 / 4 , as one would expect from 
Equation [6] 

The diurnal heating patterns shown in Figure [T] are 
similar to those measured at the surface of Earth. In 
both cases parcels of gas move in and out of the sunlight: 
on Earth this motion is entirely due to the planet's ro- 
tation (wind velocities are small compared to rotational 
velocity) , while on a gaseous planet this motion is due to 
a combination of rotation and zonal winds. Indeed, our 
model may also be relevant to a rocky planet with a thin 
atmosphere and rapid rotation, provided that the rota- 
tional period is shorter than the lateral heat conduction 
timescale. 

For a planet on a circular orbit, we can use the fact that 
the temperature extrema for a particle of gas occur when 
it reaches its equilibrium temperature. The temperature 
extrema are therefore related to their location on the 
planet by: 

f mia = cos 1 $ min (11) 

and 

Tmax = COS3 $ max , (12) 

where <& min and $ max are the angles between the sub- 
stellar meridian and temperature minimum and maxi- 

7 For simplicity, we consider here only positive e. That is, super- 
rotating equatorial jets. Later in the paper we consider both super- 
rotating and trailing winds. 
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Fig. 1. — Temperature of a parcel of gas in the atmosphere of a 
short-period planet in a circular orbit. The plot runs from dawn 
(3> = —tt/2) until dusk ($ = tt/2), and back to dawn again ($ = 
3n/2). 

mum, respectively. Since T m ; n is typically much less than 
unity, $min ~ —tt/2 in most situations (see Figure [2]), 
while T max is close to unity so small differences in this 
maximum temperature correspond to significant changes 
in the phase offset of the maximum, as shown in Figure|3l 
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-tt/2 -tt/ 4 

Fig. 2. — The solid line shows the minimum temperature reached 
by a parcel of gas, T m i n , plotted against the longitude at which 
the minimum is reached (<t> m i n = —n/2 at the dawn terminator, 
^min = at the sub-stellar point). The dashed line is the maxi- 
mum value of T m ; n = (I/71-) 1 / 4 , reached in the e — > 00 limit. We 
show the rest of the functional form of Equation [TT] as a dotted line 
for completeness. 

Note that the the diurnal heating pattern has a lon- 
gitudinal asymmetry: parcels of gas are heated faster 
than they cool, so parcels East of the hot-spot tend 
to be warmer than those West of the hot-spot. This 
asymmetry manifests itself in the disc-integrated ther- 
mal phase curve of the planet: the offset in the peak of 
the lightcurve tends to be larger than $ max . 

In the Appendix we develop analytic approximations in 
the circular regime for T max , Td us k and Td awn (the maxi- 
mum temperature reached by a parcel, its temperature at 
the dusk terminator, and at the dawn terminator). Those 
analytic approximations are compared to the numerical 
solutions in Figure [4] For intermediate recirculation ef- 
ficiencies, e « 3, the sunset temperature has a peak. 
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Fig. 3. — The solid line shows the maximum temperature reached 
by a parcel of gas, T max , plotted against the longitude at which the 
maximum is reached (<f max = at the sub-stellar point, $ maI = 
7r/2 at the dusk terminator). The dashed line is the minimum 
value of T max = (I/71-) 1 / 4 , reached in the e — » 00 limit. We show 
the rest of the functional form of Equation 1121 as a dotted line for 
completeness. 




nn r 1 1 1 1 1 

0.001 0.010 0.100 1.000 10.000 100.000 1000.000 

e 

Fig. 4. — Dimensionless temperature, T = T/To at the hot-spot, 
at dusk and at dawn, as a function of recirculation efficiency, e. 
The lines show the analytic expressions presented in the Appendix, 
while the point symbols show the result of numerically solving 
Equation [5] Although the analytic approximation does a good 
job over a variety of e, it is only applicable for planets on circular 
orbits. 

This is due to advection of heat downstream, which be- 
comes more efficient as e increases, heating the sunset 
terminator. However, when e becomes too large, heat 
is distributed uniformly along a given longitude, caus- 
ing the temperature to decrease again, hence the peak at 
intermediate e. 

4. APPLICATIONS OF MODEL 

By combining the energy transport model described 
above with a geometrical model of the planet's motions 
(both orbital and rotational), we now produce model 
phase variations for an entire planet. We model the 
planet as a grid of gas parcels: 40 longitude and 20 lat- 
itude grid points. We run the calculation for 3 orbital 
periods, with 1000 time steps per orbit. We adopt a 
perfectly edge-on orbital geometry in all cases, but oth- 
erwise use the appropriate orbital parameters (a, e, a;), 



planetary radius (i? p ), and stellar parameters (R*, T e g, 
log<7, [Fe/H]). We can run the calculations for arbitrary 
wavelengths but only present here the light curves for 
Spitzer IRAC channel 4 (8 /im). Since ours is a one- 
layer model, the shape of the phase variations is essen- 
tialy unchanged if one adopts a different waveband. The 
only free model parameters are A, r rac i, and uJ ro t] in prac- 
tice we set A — since the albedo does not significantly 
affect the shape of the thermal phase varia tions and it 
appears that hot Jupiters have low albedos ([Rowe et al.l 
20081 ICowan &: AuollEFl) 1 1). and references therein). 

All stellar and planetary data are taken from exo- 
planet.eu, maintained by Jean Schneider; using numbers 
from exoplanets.org did not perceptibly change our re- 
sults. When the stellar data are not available, we have 
assumed typical parameters for the appropriate spectral 
class, and solar metallicity. Insofar as we are only con- 
cerned with the broadband mid-IR brightnesses of the 
stars, our results should not depend sensitively on the 
input stellar parameters. Using the stars' T e g, \ogg and 
[Fe/H ], we use the PHOENIX /NextGen stellar spectrum 
grids (Hauschildt ct al. 1999) to determine their bright- 
ness temperatures at the observed frequencies. For each 
waveband, we determine the ratio of the stellar flux to 
the blackbody flux at that grid star's T e g. We then apply 
this factor to the T e g of the actual observed star. 

There are many computational shortcuts that one can 
use with this analytic model. The orbits are treated as 
edge-on and the planetary obliquity is assumed to be zero 
(this is strictly true for the core; for the atmosphere it 
simply means that the equatorial jet-stream flows in the 
East- West direction), so the D.E. need only be solved 
at one latitude on the day-side (e.g., the equator) and 
those day-side heating curves are easily adjusted for other 
latitudes via To. The D.E. can then be solved analytically 
on the planet's night-side. 

4.1. Planet- Integrated Properties 

The observed flux ratio depends on a combination of 
orbital factors (a, e), planetary factors [A, r ra( j, w r ot), 
and viewing geometry (oj and a, the usual phase an- 
gle: a = at eclipse, a = it at transit). We show 
schematically in Figure 5 how we combine the orbital, 
planetary and viewing factors to obtain disc-integrated 
thermal light curves. The top panel simply shows the 
planet's distance from its host star; the second panel 
shows the planet's orbital angular velocity; the third 
panel shows the equilibrium temperature at the sub- 
stellar point (solid line) and the highest temperature on 
the model planet (dotted line). The fourth panel shows 
the total absorbed flux (solid line) and the total emitted 
flux (dotted line). The fifth panel shows the planet's il- 
luminated fraction, / = |(1 + cos a). The bottom panel 
shows the planet/star flux ratio at 8 /im as seen from 
Earth. 

We start the calculations with all of the parcels at T = 
T , but the planet reaches a periodic equilibrium in a 
couple e- folding times (a few r ra d). If the planet had 
no heat capacity, the dotted lines would perfectly track 
the solid lines in the third and fourth panels of Figure 5. 
The effect of a non-zero planetary heat capacity (non- 
zero T ra d) i s: t ne peak-to-trough swings in emitted flux 
are damped, and the peak in emission lags the peak in 
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Fig. 5. — Diagnostics for HD 17156b, using a model with A = 0, 
T rad = 20 hours, w r ot = 2cj max . 

absorption. 

A planet on a circular orbit reaches a steady state 
whereby the temperature of a parcel of gas is always the 
same at a given latitude, 9, and local time, $. In other 
words, the planet's climate is steady as seen from the host 
star. This condition cannot occur for an eccentric planet. 
Instead, we verify that our calculations satisfy a some- 
what weaker steady-state condition: the parcels of gas 
at the same latitude, local time and orbital phase should 
have the same temperature, orbit after orbit. Since we 
use radiative timescales much shorter than any of the or- 
bital periods considered here, reaching this steady state 
is easily achieved by solving the heating curve D.E. for a 
full three orbits. 

The thermal contrast ratio between the planet and its 
host star amounts to a weighted mean of the brightness 
of the visible hemisphere of the planet: 

l=^k)(l)7 mM)B(TW (13) 

where B(Tb) is the star's intensity given its brightness 
temperature at this wavelength, B(T) is the blackbody 
intensity at the temperature T((j),9,t), and V(4>,6,t) is 
the visibility of a region on the planet (V = 1 at the sub- 
observer point, drops as the cosine of the angle from the 
observer, a nd is null on the f ar side of the planet. For 
details, see ICowan et al1l2009h . 

5. MODEL LIGHT CURVES 

Sample light curves for a planet on a circular orbit 
(HD 189733b) are shown in Figure [6] The minimum in 
the light curve occurs half an orbit from the maximum, 
despite the fact that the minimum temperature reached 
by a parcel of gas is always near the dawn terminator. 

As discussed in Agol et al. (2010), the 8 /im observa- 
tions of HD 189733b (phase variations and eclipse offset) 
indicate an e of order unity and super-rotating winds, 
shown by the blue line in Figure [6] Without further 



HD189753: P = 2.22 doys, e = 0.00 
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Fig . 6.— The 8 (im light curve for HD 189733b HBouchv et al.l 
2005) for a variety of radiative and advective timescales, and as- 
suming A = 0. The vertical dashed and dotted lines mark inferior 
conjunction (transit) and superior conjunction (eclipse), respec- 
tively. The transit and eclipses have been removed for clarity. For 
non-zero e, the min and max of the light curve will not neces- 
sarily occur at inferior and superior conjunction. The effective 
rotation rate of the atmosphere is expressed in units of the orbital 
frequency, u or b. The peak of the light curve occurs after superior 
conjunction for aj ro t < cu or h', it occurs before superior conjunction 
for OJ rot > UJ olb . 

information, it is impossible to break the degeneracy be- 
tween r r ad and w ro t (th is degeneracy has been pointed 
out in different terms in lCho e~a l. 2008a). But as men- 
tioned in § 1, many hydrodynamic al simulations of hot 
Jupiters indicate trans-sonic wi nds (jCooper fc Showman! 
120051 : iRauscher fc Menoul 120101 eg:), making for planet- 
crossing times of order 10 5 seconds (— 30 hours). This 
means that the radiative timescale of HD 189733b's at- 
mosphere is in the tens of hours. For our calculations, 
we adopt either the r ra d = limit, or a fiducial radiative 
time of 20 hours, and consider three different advective 
schemes: uj Iot = 0.5w or b is a trailing wind; it takes parcels 
of gas two planetary orbits to experience one diurnal cy- 
cle. The w ro t = 1.5cj or b is a slightly super-rotating wind; 
it again takes parcels of gas two planetary orbits to ex- 
perience one diurnal cycle. For the fast super-rotating 
winds (ui I0 t = 2.0w or b), parcels of gas experience one di- 
urnal cycle per orbit. 

The sub-rotating (w rot = 0.5cj or b) and super-rotating 
(w rot = 1.5w or b) cases have opposite phase offsets but 
identical maxima and minima. At superior conjunction 
they have identical contrast ratios, so one cannot dis- 
tinguish between the two models with secondary eclipse 
depth alone. In detail, of course, the timing and shape of 
secondary eclipse i s sensitive to the day-side brightness 
map of the planet (|Williams et al]l2006t fRausche r et al.l 
[2007tlAgol et al.1120101) . 

Figures [7HT51 show the same sample of model phase 
variations, but for planets on eccentric orbits. We use 
the same r ra d (0, 20 hrs) and u; ro t/w m ax (0-5, 1-5, 2.0) as 
for HD 189733b. For consistency and ease of comparison, 
all planets pass through periastron at t = 0, and 8 /zm 
light curves are shown. In many cases Spitzer light curves 
have been obtained at this wavelength, and since we use a 
single-layer model the shapes of the modelled lightcurves 
would be similar at other wavelengths. Furthermore, we 
assume A = across the board. Non-zero albedo would 
reduce the mean planet/star thermal contrast, but not 
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WASP-17: P = 3.74 days, e = 0.13, a = 290° 
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Fig . 7. — The 8 fim light curve for WASP-17b HAnderson et alj 
2010) for a variety of radiative and advective timescales, and as- 
suming A = 0. The vertical dashed and dotted lines mark inferior 
conjunction (transit) and superior conjunction (eclipse), respec- 
tively. The transit and eclipses have been removed for clarity. The 
effective rotation rate of the atmosphere is expressed in units of 
the maximal orbital frequency, u max . 



HAT-P-11: P = 4.89 doys, e = 0.20, a = 355° 
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Fig. 9.— The 8 £tm light curve for HAT-P-llb HBakos et al.120101) 
for a variety of radiative and advective timescales, and assuming 
A = 0. The vertical dashed and dotted lines mark inferior con- 
junction (transit) and superior conjunction (eclipse), respectively. 
The transit and eclipses have been removed for clarity. The effec- 
tive rotation rate of the atmosphere is expressed in units of the 
maximal orbital frequency, w ma x- 



GI436: P = 2.64 doys, e = 0.15, a = 551° 
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Fig. 8.— The 8 (im light curve for Gl 436b l IButler et al.ll200lh 
for a variety of radiative and advective timescales, and assuming 
A = 0. The vertical dashed and dotted lines mark inferior con- 
junction (transit) and superior conjunction (eclipse), respectively. 
The transit and eclipses have been removed for clarity. The effec- 
tive rotation rate of the atmosphere is expressed in units of the 
maximal orbital frequency, oimax- 

the shape of the light curves. Furthermore, light curves 
at multiple wavelengths are necessary to break the de- 
generacy between albedo and wavelength-dependence: a 
small DC-offset in a single waveband either indicates that 
the planet has reflected away a fraction of incident flux, 
or that the absorbed energy is being re-radiated at other 
wavelengths. 

5.1. Discussion 

The black lines (the T ra d = case) can be understood 
purely geometrically: the minimum contrast ratio of zero 
occurs at transit (marked with a vertical dashed line), 
since at that phase only the planet's night side is visi- 
ble. The maximum contrast ratio — on the other hand — 
does not necessarily coincide with superior conjunction 
(the vertical dotted line) for non-zero eccentricity. If 
the incident flux is increasing at superior conjunction 
(\oj\ < 90°), the maximum in the light curve occurs some- 



XO-3: P = 3.19 days, e = 0.26, a = 346° 
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Fig . 10.— The 8 /im light curve for XO-3b l|Johns-Krull et al.1 
2008) for a variety of radiative and advective timescales, and as- 
suming A = 0. The vertical dashed and dotted lines mark inferior 
conjunction (transit) and superior conjunction (eclipse), respec- 
tively. The transit and eclipses have been removed for clarity. The 
effective rotation rate of the atmosphere is expressed in units of 
the maximal orbital frequency, w max - 

what late (after eclipse) ; if the incident flux is decreasing 
at superior conjunction (90° < u < 270°), the peak oc- 
curs early (before superior conjunction). 

Including advection complicates this picture: super- 
rotating and trailing winds have the effect of shifting 
the peak early and late, respectively (see first Figure EJ. 
The direction of the winds at periastron (whether uj mt is 
greater than or less than w miH ) determines the offset of 
the thermal phase peak from the r ra d = case. 

We present the planets in order of increasing eccentric- 
ity, which also roughly corresponds to increasing orbital 
period. WASP- 17b (Figure[7|) has a very "circular" phase 
variation in the r ra( j = case, but the symmetry between 
the CL! ro t/w m ax = 0.5 and 1.5 cases is broken for non-zero 
radiative time. 

Gl 436b, HAT-P-llb and XO-3b (Figures M and [TOj 
show how the argument of periastron, u>, is critical in de- 
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HAT-P-2: P = 5.63 days, e = 0.52, a = 1 90 



HD80606: P = 111.78 doys, e = 0.93, a = 292° 
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Fig. 11.— The 8 fim light curve for HAT-P-2b (Bak oFetal .120071) 
for a variety of radiative and advective timescales, and assuming 
A = 0. The vertical dotted and dashed lines mark superior con- 
junction (eclipse) and inferior conjunction (transit), respectively. 
The transit and eclipses have been removed for clarity. 



HD17156: P = 21.22 days, e = 0.67, a 
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Fig . 12.— The 8 (im light curve for HD 17156b HFischer et al.l 
2007) for a variety of radiative and advective timescales, and as- 
suming A = 0. The vertical dotted and dashed lines mark supe- 
rior conjunction (eclipse) and inferior conjunction (transit), respec- 
tively. The transit and eclipses have been removed for clarity. 



termining the shape of the phase variations: the periods 
and eccentricities of the three planets differ, but they 
all have periastron shortly after superior conjunction 
and the refore strikingly simi lar li ght curves. Fo r com - 
parison, iDeming et all (|2007| ) and lDemorv et al.1 (|2007l ) 
measured 8 /xm secondary eclipse depths of 5.7(8) xlO~ 4 
and 5.4(7) xl(T 4 , respe ctively, for Gl 436b. T he multi- 
band eclipse depths of Stevenso n et al.l (|2010T l indicate 
the planet has a hot day-side; this is most consistent wit h 
the 50 x solar metalicity model of iLewis et al.1 (|2010( ). 
which exhibits a super-rotating equat orial jet. Th e 8 mi- 
cron model phase curves presented bv lLewis et al.l (|2010l ) 
are similar to our super-rotating solutions (the blue and 
green curves in Figure [8]) , and also undershoot the ob- 
served eclipse depths. 

HAT-P-2b, HD 17156b and HD 80606b (Figures Hp 
[T2l and ITS")) have such high eccentricities that even their 
Trad = lightcurves look nothing like that of a circu- 
lar planet. Furthermore, their orbital periods are suffi- 
ciently long that — for non-zero r ra d — they exhibit "ring- 



FlG. 13.— The 8 /^m light curve for HD 80606b (Na ef et al.ll200lT> 
for a variety of radiative and advective timescales, and assuming 
A = 0. The vertical dotted and dashed lines mark superior con- 
junction (eclipse) and inferior conjunction (transit), respectively. 
The transit and eclipses have been removed for clarity. Models 
run with 3x higher time-resolution yields indistinguishable light 
curves. 

ing" as the hot spot generated at periastron rotates in 
and out of view (l/u; ac iv < Trad < P)- The envelope 
is ~ exp((i — £ P eri)/7rad), an d the period of the oscil- 
lations is simply the planet's effective rotation period, 
P lot = 2ir/u! rot . In effect, the ringing caused by flash- 
heating is not dissimilar to the variability seen in some 
hydro dynamic simulations du e to precession of polar vor- 
tices (|Cho et al.l l 200ll2008al). 

For comparison. iLangton fc Laughlinl (|2008t ) predicts 
an 8 /im phase variation -Fmax/^min of 2.8 for HAT- 
P-2b, similar to our values of ~ 2 for super-rotating 
winds (see T able [T] for a summ ary of our results). For 
HD 1715 6b. Iirwin et al l (120081) predicts F max /F min of 
2.9, while llrofc Dcmingi (|2010| ) predicts 3.3; our models 
predict somewhat greater phase variations, .F m ax/-Fmin ~ 
5. Note that if periastron occurs near inferior conjunc- 
tion, our model predicts that the maximum contrast ra- 
tio may actually be greater for non-zero T ra d (e.g., Fig- 
ure [H]). 

Observations of HD 80606b show an 8 ^m eclipse depth 
of 1.0(2) xl0~ 3 (|Laughlin et al.ll20 0i). similar to our sub- 
rotati ng model (red line i n Fig ure I13|) . For compar- 
ison, Langton & Laughlin (2008) predicts phase varia- 
tions with -F m ax/-Fmin = 2.1, much smaller that our val- 
ues of ~ 60. Note that they have r ra d = 4.5(2) hours 
at the 8 /xm photosphere, but the planet never cools to 
below a contrast ratio of ~ 4 x 10~ 4 , despite the planet's 
112 day orbit. This discrepancy is largely due to our im- 
plicit assumption of no internal heat: after spending over 
100 days far from its host star, our model planet cools 
to ~ 400 K. This is still much hotter than the remnant 
heat of formation (50-100 K), but cooler than the 720 K 
night-side equilibriu m temperature used for the h ydrody- 
namic simulations offLaneton fe Laughlinl (|2008l ) . Mean- 
while, llro fc Demine) (|2010f ) predicts F m ,v/F mi „ = 4.2. 
The peak to trough a mplitudes of our mode ls are w ithin 
a factor of 2 of th e ILangton fc Laughlinl (|2008f ) and 
llro fc Deminel (120101) simulations. If we ado pt the r ra d = 
4.5(2) hour of ILangton fc Laughiml (|2008f ). we obtain 
peak-to-trough amplitudes of ~ 3 x 10~ 3 , still a factor of 
~ 2 larger than their models. That factor of 2 is signifi- 
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cant: it indi cates that — unlike th e lLangton fe Lau ghlin 
(|2008ft and llro fc Demingl (|2010fl simulations— m our 
model the planet can completely radiate away it's heat 
over a single planetary orbit. 

In general, non-zero heat capacity reduces the max- 
imum contrast ratio and increases the minimum, mak- 
ing for a flatter light curve. Furthermore, the phase 
amplitude generally diminishes for larger advective fre- 
quencies: the peak-to-trough amplitude is smaller for 
w ro t = l-5w max than for Lu Iot = 2w max across the board. 
For planets with 90° < cj < 270°, however, a second or- 
der effect may also take place: the light curve maximum 
increases for moderate T ra d as this allows the planet to 
hold onto absorbed power until its hot-spot has rotated 
into our view (see Figure [T2|) . 

6. MODEL ASSUMPTIONS 
Agnostic about Tidal Locking: Short-period exoplan- 
ets on circular orbits should have tidally locked cores, 
because the tidal locking timescale is considerably 
shorter than the age of most planetary systems (e.g., 
ILubow et aill 997). By the same token, planets on eccen- 
tric orbits should have cores in pseudo-synchronous ro- 
tation (roughly speaking, this means being tidally locked 
at periastron, when the tidal forces are strongest). Un- 
like detailed hydro simulations that use the planet's ro- 
tation rate as an input parameter, our model is agnostic 
about precise prescripti ons for pseudo- synchronous rota- 
tion frequency, w ps (e.g. lHutlll981t Hvanov fc Papaloizoul 
120071) . For e x ample , at a moderate eccentricity of e = 
0.3, the iHutl (1 19811) formulation give s w ps « 0.8w max , 
while the llvanov fc Papaloizoul (|2007f) formulation gives 
cu ps w 1.4w max . These details about the planet's interior 
are not important for our purposes, however, because the 
planets in our study have thick gaseous envelopes cover- 
ing their cores and hydrodynamic simulations indicate 
that wind-speeds in the vicinity of the photosphere are 
comparable to the rotational velocity. As a result, the 
parcels of gas are probably not stationary with respect to 
the sub-stellar point, even for (presumably tidally locked) 
planets on circular orbits. As a first-order model we treat 
them as advecting longitudinally at a constant angular 
velocity. 

No Internal Heating: Tidal circularization is undoubt- 
edly responsible for the generally low ecce ntricity of 
short-period planets (e.g. iFord fc Ra sio 2008, and refer- 
ences therein). It is not clear, however, that tidal heating 
has an observable e ffect on light curves. For example, 
Uackson et al.l (|2008l ) calculate present-day tidal heating 
between 2 and 5 x 10 17 W for Gl 436b, two orders of 
magnitude less power than the stellar insolation, which 
ranges between 7 and 13 x 10 19 W. For short-period plan- 
ets, the observational signature of internal heat (either 
remnant gravitational energy from formation or ongoing 
tidal heating) will likely be an inflated radius rather than 
additional emergent flux. If hot Jupiters exist for which 
the tidal heating is a sizeable fraction of the insolation, 
then our model will have under-estimated the zeroeth- 
order (DC) component of their lightcurves. By the same 
token, we also neglect the remnant heat of formation, 
which is likely less t han 100 K for a Gyrs-old planet 
(|Burrows et al.l I2006D . One important consequence of 
neglecting any internal heating and having a simple one- 
layer model, is that the night-side equilibrium temper- 



ature is zero. For the most part this is not important, 
since parcels pass through the illuminated hemisphere 
often enough to never get very cool. The exception is 
for our models of HD 80606b, where the planet spends 
such a long time (> 100 days) far from its host star. It 
is possible that the real planet never cools to the 400 K 
seen in our toy model. 

No Heat Conduction: Our model neglects conduction 
of heat from one cell to its neighbors. As such, our model 
is only relevant when advection and radiation are the 
dominant energy transport mechanisms. This is proba- 
bly the case for hot Jupiters, which are thought to have 
strong winds. If heat conduction is important on hot 
Jupiters, then they may exhibit smaller day-night tem- 
perature contrasts and the relation between the thermal 
phase offset and the thermal phase amplitude will be 
broken (see first Equation 12). 

No Mechanical Work: Our model assumes that all the 
incident energy from the host star is either reflected or 
absorbed, and the energy that is absorbed only goes 
into increasing the temperature of the parcels of gas. 
In effect, this amounts to neglecting the work term in 
the First Law of Thermodynamics. In the steady-state, 
the power absorbed by the planet must be balanced by 
the power emitted. That is to say, the ultimate ef- 
fect of work is to warm up parcels of gas (e.g., parcels 
speed up to su per-sonic velociti es, shock and deposit 
energy as heat; IGoodm an 2009). This implicitly ne- 
glects four effects that may very well be occurring: en- 
ergy may instead go i nto latent heat to vapor ize partic- 
ulates; thermal tides (jArras fc Socrates! 120101); increas- 
ing t he atmospheric scale height ()Guillot fc Showman! 
2002); or speeding up winds. Any of these "work" 
terms would tend to dampen the temperature fluctua- 
tions experienced by a parcel of gas: they act as en- 
ergy sinks near periastron and may act as energy sources 
near apoastron. As such our model will under-estimate 
gas response time, over-estimate the temperature excur- 
sions, and hence the peak-to-trough amplitude of disc- 
integrated lightcurves. The planet most likely to show 
these effects is HD 80606b, since realistically much of the 
energy absorbed at periastron must go into speeding up 
winds, inflating the atmosphere, etc. Indeed, the peak- 
to-trough phase amplitudes we compute for this planet 
are roughly a factor o f 2 greater than those from hydro- 
dynamic simulations (Langton & Laughlin 20081) . 

No Latitudinal Advection: Our toy model neglects ad- 
vection towards the planet's poles: effectively we have 
zonal but no meridional winds. At low (mbar) pressures, 
there almost certainly is meridional flow away from the 
sub-stellar point, but the low presures mean that only a 
small fraction of absorbed power is moved in this way. 
But even when meridional flows are not efficient, lat- 
itudinal heat tran sport may still o c cur in the form of 
eddy fluxes (e.g., ICho et all [20081 IChol l2008h . which 
we have also implicitly neglected. Efficient equator-pole 
heat transport on a planet would have the effect of re- 
ducing the DC component of the phase curve. In the 
limit of perfect latitudinal heat transport, the equatorial 
temperature is suppressed by a factor (tt/4) 1 / 4 w 0.94 
and the disc-integrated temperature is only suppressed 
by (37r 2 /32) 1 / 4 w 0.98 compared to the nominal (no 
latitudinal transport) case. This latitudinal dilution is 
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(necessarily) smaller than the longitudinal dilution. As 
pointed out in lBudail (|2010f ). latitudinal energy transport 
is only important for exoplanets that are not edge-on (we 
get a better view of their pole), and in those cases the 
unknown planetary radius is a confounding factor. 

Single Wind Velocity: Our model uses a single angular 
wind velocity, which can roughly be thought of as the ve- 
locity of the equatorial jet stream. We are neglecting the 
jet-streams at mid-l atitudes that transport en ergy in the 
opposite direction (Dobbs- Dixon et al.ll2010f) . We have 
assumed that — to first order — the velocity of the equa- 
torial jet is the only speed that affects the observed phase 
variations of a p lanet, but the speed of waves may also be 
important (e.g. JWatkins fc Cho |[2010h . In any case, our 
w ro t can be thought of as the visibility-weighted mean of 
the competing zonal jets, which is necessarily slower than 
the velocity of the equatorial jet. Note that our adoption 
of constant an gular velocity a mounts to solid body ro- 
tation (see also llro et al.M2005h . As far as simplifications 
go, solid-body rotation is attractive because it naturally 
produces smaller linear velocities near the planet's poles 
than near the equator, a generic feature of many hydro- 
dynamic simulations. That being said, the assumption 
is undoubtedly wrong in detail and one could instead 
parameterize the zonal flow as having a constant linear 
rather than angular velocity: w rot (6 l ) = w ro t sin 9. Such a 
parameterization could be normalized to yield the same 
overall heat transport, but it would smear out thermal 
phase curve features, since parcels of gas flash- heated at 
periastron would not move in lockstep across the face of 
the planet. 

Constant Wind Velocity: We assume that the wind 
speed does not change with time. The existence 
of steady-state, E-W jet streams is supported by 
hydrody namic simulatio n s of pl anets on circular 
orbits dShowman et al.l l2009at iRauscher fc Menoul 
I2010t iDobbs-Dixon fc Lin! 12*0081 and re fe rences 
there in). But the simu l ations of iCho et all (2003) 
and pThrastar son fc Cho I (|2010f ) indicate that these 
results may be a function of initial conditions, in 
which case even planets on circular orbits may exhibit 
wandering equatorial jets. For planets on eccentric 
orbits, the wind speed of the equatorial jet may very 
well slow down as the planet's energy budget decreases, 
but the timing of the primary peak should only depend 
on the wind speed near pericenter. If the winds damp 
out near apastron, then the maximum wind velocity 
may occur shortly after periastron, depending on the 
inertia of planet's atmosphere. In effect, our model 
works in the limit of large inertia, so that the winds do 
not speed up or slow down even as the planet's energy 
budget changes. There is very little energy to transport 
when the planet reaches apoastron, so the wind speed 
at that point in the orbit does not have observable 
consequences. 

Constant Albedo: Atmospheric albedo is thought to in- 
crease with decreasing temperature. This can be quite 
sudden if, for example, clouds form high in the atmo- 
sphere (eg. lKane fc Gelinoll2010l and reference s therein). 
As mentioned in iLangton fc Laughlinl (|2008bl V the for- 
mation of clouds on an eccentric planet as it moves away 
from periastron will simply exacerbate its 1/r 2 energy 
budget, making the periastron passage even more im- 



portant for thermal light curves. Note that this is the 
only of our model assumptions that could lead us to un- 
rfer-estimate the amplitude of phase variations. 

Blackbody Emitters: While we treat the parcels of gas 
as blackbody emitters, the shapes of our lightcurves 
should not change much if the planet has a non- 
blackbody spectrum. The shape of the light curves is 
preserved provided that the mid-IR flux increases as a 
simple function of temperature. This is likely to be true 
on the Rayleigh- Jeans tail of the planet's spectral energy 
distribution. 

Vertically Isothermal: Using a single temperature at 
each position on the planet (ie: a single-layer model) 
implicitly means assuming that the temperature at the 
emitting surface is roughly equal at all wavelengths of 
interest. In the case of highly irradiated planets, the 
day-side temperature-pressure profile is shallower than 
for an isolated giant planet. Depending on how much 
advection of heat occurs — and at what depths — this 
may not be true of the night-side. Observationally, the 
day-side brightness temperatures of hot Jupiters often 
differ from one infrared waveband to another. But the 
only hot Jupiter with thermal phase curves observed at 
multiple wavelengths is HD 189733b, and the shape of 
the 8 and 24 micron ph ase curves are indistin guishable, 
given the uncertainties (jKnutson et al.l l2009b|) . This is 
notable since 1-D atmospheric models of this planet place 
the 8 micron photosphere at more than twice the pres- 
sure of the 24 micron photosphere. Note that this co- 
incidence for HD 189733b only means that the ratio of 
the radiative and advective times is the similar at the 8 
and 24 micron photospheres. For eccentric planets the 
degeneracy between those two timescales is broken, so 
seeing the same morphology at two different wavebands 
means that both the radiative and advective timescales 
are the same at the two photospheres. If the phase vari- 
ation morphology is grossly different from one waveband 
to another, then one would need a model with different 
radiative and advective times as a function of waveband 
(or equivalently, a multi-layer model). 

7. CONCLUSIONS 

We have presented a semi-analytic model for the heat- 
ing pattern of gas in the atmosphere of a hot Jupiter. Our 
two physically-motivated parameters are the planet's ra- 
diative timescale, and its wind speed. This model is 
meant to calculate the shape of thermal phase variations 
of transiting exoplanets with any eccentricity or argu- 
ment of periastron. We have assumed zero albedo and 
blackbody- like emission throughout: these assumptions 
have an effect on the DC offset of our phase variations 
(e.g., the secondary eclipse depth), but should not signif- 
icantly affect the shape of the phase variations (timing 
of the phase peak, peak-to-trough amplitude, etc.). 

Our model is sufficient to fit the shape of phase vari- 
ations for HD 189733b, a hot Jupiter on a circular or- 
bit. The radiative and advective times of our "best fit" 
model for that planet are within the range of plausible 
parameters from detailed hydrodynamic simulations. We 
bracket these parameters and compute light curves for a 
selection of eccentric transiting planets. Our calculated 
8 /zm light curves are summarized in Table [TJ 

Our principal findings from this numerical experiment 
are: 
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TABLE 1 

Calculated Phase Variations of Eccentric Planets 



r^lanet 




ad 


<^rot / "max 


Amplitude 6 


At c 

max 


HD 18Q7SSb 

i 1 1—> i O ,J J lj l_f 


o 


hrs 


N/A 


4.9 


x 


10~ 3 


0.0 


hrs 




20 


hrs 


0.5 






2.9 


5.3 


hrs 




20 


lll'S 


1.5 






2.9 


-5.3 


hrs 




20 


hrs 


2.0 






1.9 


-7.2 


hrs 


WASP-17b 





Ins 


N/A 


4.9 


x 


10 -3 


0.27 


hrs 




20 


hrs 


0.5 






3.9 


5.8 


hrs 




20 


hrs 


1.5 






2.4 


-2.4 


hrs 




20 


lll'S 


2.0 






1.8 


-4.7 


hrs 


d] 436h 


o 




N/A 


8.0 


x 


10~ 4 


1.8 


hrs 




20 


hrs 


0.5 






6.5 


5.0 


hrs 




20 


lll'S 


1.5 






2.4 


-24 


hrs 




20 


hrs 


2.0 






1.7 


-31 


hrs 


HAT-P-llb 





hrs 


N/A 


5.1 


X 


10~ 4 


2.6 


hrs 




20 


hrs 


0.5 






9.4 


6.7 


hrs 




20 


hrs 


1.5 






2.7 


-2.2 


hrs 




20 


hrs 


2.0 






1.8 


-4.7 


lll'S 


XO-3b 





hrs 


N/A 


2.6 


x 


10 -3 


1.6 


hrs 




20 


hrs 


0.5 






3.9 


3.8 


hrs 




20 


hrs 


1.5 






1.5 


-42 


hrs 




20 


hrs 


2.0 






1.4 


-18 


hrs 


HAT-P-2b 







N/A 


2.1 


X 


10~ 3 


-3.9 


hrs 




20 


hrs 


0.5 






4.1 


3.0 


hrs 




20 


hrs 


1.5 






2.4 


-5.1 


hrs 




20 


hrs 


2.0 






2.1 


-6.6 


hrs 


HD 17156b 





hrs 


N/A 


5.9 


X 


10- 4 


-82 


hrs 




20 


hrs 


0.5 






4.9 


-68 


hrs 




20 


hrs 


1.5 






5.4 


-95 


hrs 




20 


hrs 


2.0 






5.3 


-99 


hrs 


HD 80606b 





hrs 


N/A 


3.0 


X 


10~ 3 


1.8 


hrs 




20 


hrs 


0.5 






63.5 


7.2 


hrs 




20 


hrs 


1.5 






56.1 


26 


hrs 




20 


hrs 


2.0 






51.0 


21 


hrs 



5. Planets with non-zero radiative times have their 
thermal phase peaks offset early for super-rotating 
winds, and late for sub-rotating winds. 

6. For planets with non-zero r ra d, the peak-to-trough 
phase amplitude generally decreases with increas- 
ing w ro t; for planets on circular orbits this can be 
understood in terms of the dimensionless advective 
efficiency, e. 

7. For planets that are approaching their star at su- 
perior conjunction or planets on circular orbits, in- 
creasing r ra( j has the effect of reducing the light 
curve maximum. 

8. Planets that are approaching their star at supe- 
rior conjunction also exhibit stronger phase varia- 
tions if they have sub-rotating (rather than super- 
rotating) winds at periastron. This is because the 
hot spot created at periastron rotates out of view 
more slowly. 

9. For planets receding from their star at superior con- 
junction, an increased r rac i may increase the light 
curve maximum. 

10. Eccentric planets with orbital periods significantly 
longer than T ra d exhibit "ringing" whereby the hot 
spot generated at periastron rotates in and out of 
view. 



a Wc a dopt A — for a ll these c a lcul a tions. The sh a pe of the 
phase v a ria,tions is not a ffcctcd by a lbcdo, to first order. 

For T ra d — 0, this refers to the oiaximum flux r a tio, Fmax/F*: 
the minimum flux ratio is always zero in th a t c a se, so the peak- 
to-trough a mplitudc is cqu a l to the peak v a luc. For the non- 
zero r ra d, the amplitude refers to the ratio of the maximum and 
minimum of the light curve, i^max/^min- 

c Time offset between the maximum planet/star contrast and 
superior conjunction. Negative values mean that the maximum 
occurs before eclipse; positive values mean that maximum occurs 
after eclipse. 

1. For planets on circular orbits, the location on the 
planet of the primary hot spot is intimately related 
to its amplitude. 

2. The minimum temperature of a parcel occurs at or 
shortly after dawn, but the minimum in the disc- 
integrated phase variations occurs when the pri- 
mary hotspot is on the far side of the planet rather 
than when the dawn terminator is facing the ob- 
server. 

3. For a circular orbit, the timing of the phase varia- 
tion maximum with respect to superior conjunction 
indicates the direction of the dominant winds, but 
cannot break the degeneracy between wind speed 
and radiative time. 

4. For eccentric planets — in the limit of no 
advection — the light curve maximum leads or 
trails superior conjunction depending on whether 
the planet is receding from or approaching the star 
at superior conjunction. 



11. The existence of ringing makes it possible to ex- 
plicitely measure both w rot (the frequency of the 
ringing) and T ra d (the damping of the ringing). 

As discussed in § 6, the myriad assumptions made in our 
model will tend to exaggerate the amplitude of thermal 
phase variations. This could lead, for example, to over- 
estimating the radiative time if one uses our models to 
interpret observed phase variations. 

An obvious direction for the future is to couple our 
simple the rmodynamic model to a 1-D radiative transfer 
code (e.g. Ilro fc De ming 20l3) to produce time-variable 
spectra of eccentric planets rather than just light curves. 
Such a hybrid solution would require far fewer computa- 
tional resources than full time- variable radiative transfer. 
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APPENDIX 

In the slow cooling regime (e 3> 1), the changes in 
temperature during each orbit are small, so one can 
treat Equation 10 perturbatively. The trial solution is 
T = T + ST, where T = 7r _1 / 4 is a constant equal to 
the average temperature and ST contains all of the vari- 
ability. Substituting this into Equation [10] and noting 
that max(cos $, 0) = |(cos<I> + |cos<E>|), we get: 



dT l f l < * 
— = - - cos $ 

d$ e V 2 1 



cos<E>|] - - 

7T 



AT^ST ) , (14) 



where we have only kept terms up to first-order in ST. 
The solution to this equation is 



~ _ 3^ | 7cos$ + sm$ 
day "4 ° + e(l + 7 2 ) 



-7* 



knight — jT ■ 



-7($ — 7r) 



2e(l + 7 2 )sinh (777/2) 



2e(l +7 2 )sinh (777/2)' 



Note that the day-side phase peak (middle term) has 
the right limits: for large e, the variation is as sin $ (slow- 
cooling regime), while for small e, the variations is as 
cos<& (instant-cooling limit). 

In the circular case, the longitude of the peak of the 
temperature, $ max , goes from for e = to 7r/2 for 
e — > 00 (but note that in practice the parcel's heating 
curve becomes flat before the the peak reaches tt/2, as 
shown in Figure [3]). We compute $ max numerically, and 
fit tan($ max (e)) with the function 



f(e) = x (l + Xl e- X2+{xi/{1+Xi,i)) )- 1 1 



for which we find a good fit for xq = 2.9685, x\ = 7.0623, 
x 2 = 1.1756, x 3 = -0.2958, and x 4 = 0.1846. This then 
gives the maximum temperature as a function of e as 

T max (e) = cos 1 / 4 $ max (e) « cos 1 / 4 (tan" 1 /(e)) . (17) 

The maximum temperature goes from 1 at small e (in- 



stant re-radiation) to tt" 1 / 4 at large e (uniform tempera- 
ture) . We have computed T max numerically, shown with 
the asterixes in Figure [4] while Equation [17] is plotted 
as the solid line; as can be seen the agreement is quite 
good. 

The next critical characteristics are the temperatures 
at dawn and sunset — these determine the temperature 
of the atmosphere at the terminator, where it absorbs 
starlight during transit. We can estimate the sunset tem- 
perature as follows. For very small values of e, the flux 
absorbed on the day-side is emitted nearly instantly, so 
heating balances cooling, and Td ay (e «C 1) ~ cos 1 / 4 *! 1 . 
However, near sunset, $ = 7r/2, the temperature drop 
becomes so significant that dT day /d$ becomes compa- 
rable to the cooling rate, —008$. At this point the 
planet's temperature is no longer determined by instant 
re-radiation of the absorbed heat; instead we need to take 
into account the cooling timescale. For small e this occurs 
very close to tt/2, so we can estimate the temperature at 
sunset as the time when edT day (e < 1) « — cos$. This 
(15) gives edfeos 1 / 4 $]/d$ = - \e cos~ 3 / 4 $ sin $ « -cos<I>. 

Since $ w 7r/2 and f « cos 1 / 4 $, this gives Tdusk ~ 
(e/4) 1 / 7 . This is a very weak dependence on e. For 
large values of e, recirculation is very efficient, so Td us k ~ 
^-1/4 Combining these two limits, we find the approx- 
imate formula f dusk « [tt 2 (1 + y /e)- s + y^ 8 / 7 ]" 1 / 8 , 
where yo = 0.69073 and y\ = 7.5534, which is plotted as 
the dotted line in Figure 2] as well as numerically com- 
puted data points (diamonds); again the agreement is 
quite good. 

Finally, the minimum in temperature occurs at (or very 
close to) the dawn terminator after the gas has cooled 
all night and before it starts warming again on the day 
side. The dawn temperature is well approximated by 
Tdawn = [tt + (3n/e) 4/3 }- 1/4 , plotted as the dashed line 
in Figure 2] as well as numerically computed data points 
(triangles). The agreement for this expression is quite 
accurate as well. 



(16) 
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